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Abstract — A new analysis is presented for the direct-sequence 
code-division multiple access (DS-CDMA) cellular uplink. For a 
given network topology, closed-form expressions are found for 
the outage probability and rate of each uplink in the presence 
of path-dependent Nakagami fading and log-normal shadowing. 
The topology may be arbitrary or modeled by a random spatial 
distribution for a fixed number of base stations and mobiles 
placed over a finite area with the separations among them 
constrained to exceed a minimum distance. The analysis is 
more detailed and accurate than existing ones and facilitates 
the resolution of network design issues, including the influence 
of the minimum base-station separation, the role of the spreading 
factor, and the impact of various power-control and rate-control 
policies. It is shown that once power control is established, the 
rate can be allocated according to a fixed-rate or variable-rate 
policy with the objective of either meeting an outage constraint 
or maximizing throughput. An advantage of the variable-rate 
policy is that it allows an outage constraint to be enforced on 
every uplink, whereas the fixed-rate policy can only meet an 
average outage constraint. 

I. Introduction 

The classical analysis of the uplinks (e.g., Q], (2), 0) in a 
cellular network entails a number of questionable assumptions. 
Among the principal ones are the existence of a lattice or 
regular grid of base stations and the modeling of intercell 
interference at a base station as a fixed fraction of the total in- 
terference. Although conceptually simple and locally tractable, 
the grid assumption is a poor model for actual base-station 
deployments, which cannot assume a regular grid structure 
due to a variety of regulatory and physical constraints. The 
intercell-interference assumption is untenable because the frac- 
tional proportion of intercell interference varies substantially 
with the mobile and base-station locations, the shadowing, and 
the fading. 

More recent analyses of cellular networks (e.g., ||4), |]5]) 
locate the mobiles and/or base stations according to a two- 
dimensional Poisson point process (PPP) over a network that 
extends infinitely on the Euclidian plane, thereby allowing 
the use of analytical tools from stochastic geometry [6|, 
Q. Although the two principal limitations of the classical 
approach are eliminated, the PPP approach is still unrealistic 
because it cannot model networks of finite area and does 
not permit a minimum separation between base stations or 
between mobiles, although both minimum separations are 
characteristic of actual macro-cellular deployments. 
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The analysis in this paper is driven by a new closed-form 
expression J8) for the conditional outage probability of a 
communication link, where the conditioning is with respect to 
an arbitrary realization of the network topology. The appro ach 
involves drawing realizations of the network according to the 
desired spatial and shadowing models, and then computing the 
outage probability. Because the outage probability is averaged 
over the fading, it can be found in closed form with no need 
to simulate the corresponding channels. A Nakagami-m fading 
model is assumed, which models a wide class of channels, 
and the fading parameters do not need to be identical for 
all communication links. This flexibility allows the modeling 
of distance-dependent fading, where mobiles close to the 
base station have a dominant line-of-sight (LOS) path, while 
the more distant mobiles are non-LOS. The analysis can 
characterize more than just the average performance of a 
typical uplink. The outage probability of each uplink can be 
constrained, and the statistics of the rate provided to each 
mobile for its uplink can be determined under various power- 
control and rate-control policies. 

Although our analysis is applicable for arbitrary topolo- 
gies, we focus on a constrained random spatial approach to 
model the direct-sequence code-division multiple access (DS- 
CDMA) uplink. The spatial model places a fixed number 
of base stations within a region of finite extent. The model 
enforces a minimum separation among the base stations for 
each network realization, which comprises a base-station 
placement, a mobile placement, and a shadowing realization. 
The model for both mobile and base-station placement is a 
binomial point process (BPP) with repulsion, which we call a 
uniform clustering model [8|. 

A DS-CDMA uplink differs from a downlink, which has 
been presented in a separate paper [9], in at least four 
significant ways. First, the sources of interference are many 
mobiles for the uplink whereas the sources are a few base 
stations for the downlink. Second, sectorization is a critical 
factor in uplink performance whereas it is of minor importance 
in downlink performance. Third, the uplink signals arriving at 
a base station are asynchronous and, hence, non-orthogonal. 
Fourth, base stations are equipped with better transmit high- 
power amplifiers and receive low-noise amplifiers than the 
mobiles. The net result is that the operational signal-to-noise 
ratio is typically 5 - 10 dB lower for the uplink than for the 
downlink. 
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Fig. 1. Example network topologies. Base stations are represented by large 
red circles, and cell boundaries are represented by solid black lines. Left sub- 
figure: Actual base station locations from a current cellular deployment. Right 
subfigure: Simulated base-station locations using a base-station exclusion zone 
r bs = 0.25. In the right subfigure, the simulated positions of the mobiles are 
represented by black dots, the sector boundaries are represented by light blue 
lines, and the average cell load is M/C = 16 mobiles. 

II. Network Model 

Cells may be divided into sectors by using several di- 
rectional sector antennas or arrays, each covering disjoint 
angular sectors, at the base stations. Only mobiles in the di- 
rections covered by a sector antenna can cause multiple-access 
interference on the uplink from a mobile to its associated 
sector antenna. Thus, the number of interfering signals on 
an uplink is reduced approximately by a factor equal to the 
number of sectors. Practical sector antennas have patterns with 
sidelobes that extend into adjacent sectors, but the performance 
degradation due to overlapping sectors is significant only for 
a small percentage of mobile locations. Three ideal sector 
antennas and sectors per base station, each covering 27r/3 
radians, are assumed in the subsequent analysis. The mobile 
antennas are assumed to be omnidirectional. 

The network comprises C base stations and cells, 3C 
sectors {Si, ..., S^c}, and M mobiles {X\, ..., Xm}- The base 
stations and mobiles are confined to a finite area, which is 
assumed to be a circle of radius r net and area Ttr% et . The 
sector boundary angles are the same for all base stations. 
The variable Sj represents both the j th sector antenna and 
its location, and JQ represents the i th mobile and its location. 
An exclusion zone of radius rbs surrounds each base station, 
and no other base stations or mobiles are allowed within this 
zone. Similarly, an exclusion zone of radius r m surrounds each 
mobile where no other mobiles are allowed. 

Example network topologies are shown in Fig. Q] The left 
subfigure shows the locations of actual base stations in a small 
city with hilly terrain. The base-station locations are given by 
the large filled circles, and the Voronoi cells are indicated in 
the figure. The minimum base-station separation is observed 
to be about 0.43 km. The right subfigure shows a portion of a 
randomly generated network with average number of mobiles 
per cell M/C — 16, a base-station exclusion radius r^ s = 
0.25, and a mobile exclusion radius r m = 0.01. The locations 
of the mobiles are represented by small dots, and light lines 
indicate the angular coverage of sector antennas. 

Consider a reference receiver of a sector antenna that 
receives a desired signal from a reference mobile within its 
cell and sector. Both intracell and intercell interference are 
received from other mobiles within the covered angle of the 



sector, but interference from mobiles in extraneous sectors is 
negligible. Duplexing is assumed to prevent interference from 
other sector antennas. The varying propagation delays from 
the interfering mobiles cause their interference signals to be 
asynchronous with respect to the desired signal. 

In a DS-CDMA network of asynchronous quadriphase 
direct-sequence systems, a multiple-access interference signal 
with power / before despreading is reduced after despreading 
to the power level Ih(r )/G, where G is the processing gain or 
spreading factor, and h(r ) is a function of the chip waveform 
and the timing offset r of the interference spreading sequence 
relative to that of the desired or reference signal. If r is 
assumed to have a uniform distribution over [0, T c ], then the 
expected value of h(r ) is the chip factor h . For rectangular 
chip waveforms, h = 2/3 ifTOl . ifTTI . It is assumed henceforth 
that G/h(r ) is a constant equal to G/h at each sector receiver 
in the network. 

Let Aj denote the set of mobiles covered by sector antenna 
Sj. A mobile Xi G Aj will be associated with Sj if the 
mobile's signal is received at Sj with a higher average power 
than at any other sector antenna in the network. Let Xj C 
Aj denote the set of mobiles associated with sector antenna 
Sj. Let X r 6 Xj denote a reference mobile that transmits a 
desired signal to Sj. The power of X r received at Sj is not 
significantly affected by the spreading factor and depends on 
the fading and path-loss models. The power of Xj, i ^ r, 
received at Sj is nonzero only if Xi G Aj, is reduced by 
G/h, and also depends on the fading and path-loss models. 
We assume that path loss has a power-law dependence on 
distance and is perturbed by shadowing. When accounting for 
fading and path loss, the despread instantaneous power of Xi 
received at Sj is 
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(1) 
where gij is the power gain due to fading, £ij is a shadowing 
factor, Pi is the power transmitted by Xi, and Aj\X r is set 
Aj with element X r removed (required since X r does not in- 
terfere with itself). The {gij} are independent with unit-mean, 
and g 
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where a^j is Nakagami with parameter rriij. 
While the {gij} are independent from each mobile to each 
base station, they are not necessarily identically distributed, 
and each link can have a distinct Nakagami parameter rriij. 
When the channel between Sj and Xi experiences Rayleigh 
fading, m^j = 1 and gij is exponentially distributed. In the 
presence of log-normal shadowing, the {£,ij} are i.i.d. zero- 
mean Gaussian random variables with variance a\. In the 
absence of shadowing, £j j = 0. The path-loss function is 
expressed as the attenuation power law 
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where a > 2 is the attenuation power-law exponent, and do 
is sufficiently large that the signals are in the far field. It is 



assumed that r^ s and r m exceed do- 

It is assumed that the {g%,j} remain fixed for the duration of 
a time interval, but vary independently from interval to interval 
(block fading). With probability pi, the i th interferer transmits 
in the same time interval as the reference signal. The activity 
probability {pi} can be used to model voice-activity factors or 
controlled silence. Although the {pi} need not be the same, it 
is assumed that they are identical in the subsequent examples. 

Let g(i) denote a function that returns the index of the sector 
antenna serving X{ so that Xi £ Xj if g (i) = j. Usually, the 
sector antenna S s u-\ that serves mobile Xi is selected to be 
the one with index 

g (?) = argmax |l0«-/ 10 / (||5, - X t \\) , X t G 

which is the sector antenna with minimum path loss from Xi 
among those that cover Xi. In the absence of shadowing, it 
will be the sector antenna that is closest to Xi. In the presence 
of shadowing, a mobile may actually be associated with a 
sector antenna that is more distant than the closest one if the 
shadowing conditions are sufficiently better. 

The instantaneous SINR at sector antenna Sj when the 
desired signal is from X r £ Xj is 

Pr,j 
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where J\f is the noise power and Ii is is a Bernoulli variable 
with probability P[I t = I] = pi and P[Ii = 0] = 1 - p % . 
Substituting {Q) and (fJJ) into © yields 
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where T = d^P r /N is the signal-to-noise ratio (SNR) due to 
a mobile located at unit distance when fading and shadowing 
are absent, and 
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is the normalized mean despread power of Xi received at 
Sj, where the normalization is by P r . The set of {&i,j} 
for reference receiver Sj is represented by the vector ftj = 
{Oij, ...,Qmj}- 

III. Outage Probability 

Let p r denote the minimum instantaneous SINR required 
for reliable reception of a signal from X r at its serving sector 
antenna Sj,j — g(r). An outage occurs when the SINR of a 
signal from X r falls below (3 r . The value of j3 r is a function of 
the rate R r of the uplink, which is expressed in units of bits per 
channel use (bpcu). The relationship R r = C(/3 r ) depends on 



the modulation and coding schemes used, and typically only a 
discrete set of R r can be selected. While the exact dependence 
of R r on f3 r can be determined empirically through tests 
or simulation, we make the simplifying assumption when 
computing our numerical results that C((3 r ) — log 2 (l + {3 r ) 
corresponding to the Shannon capacity for complex discrete- 
time AWGN channels. This assumption is fairly accurate for 
modern cellular systems such as Wideband CDMA (WCDMA) 
or LTE, which use turbo codes with a large number of available 
rates. 

Conditioning on ilj, the outage probability of a desired 
signal from X r £ Xj that arrives at Sj is 



e r = P ry r ,j < Pr Oj] 



(7) 



Because it is conditioned on ilj, the outage probability 
depends on the particular network realization, which has 
dynamics over timescales that are much slower than the fading. 
In (8), it is proved that 
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IV. Network Policies 



A. Power Control 



A standard power-allocation policy for DS-CDMA networks 
is to select the transmit power {Pi} for all mobiles in the set 
Xj such that, after compensation for shadowing and power-law 
attenuation, each mobile's transmission is received at sector 
antenna Sj with the same power Pq. For such a power-control 
policy, each mobile in Xj will transmit with a power Pi that 
satisfies 



P i ltf<j' 1Q f(\\S j -X i \\) = P , x % 
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Since the reference mobile X r £ Xj, its transmit power is 
determined by (fT2l . 

For a reference mobile X r , the interference at sector antenna 
Sj is from the mobiles in the set Aj\X r . This set can be 
partitioned into two subsets. The first subset Xj\X r comprises 
the intracell interferers, which are the other mobiles in the 
same cell and sector as the reference mobile. The second 
subset Aj\Xj comprises the intercell interferers, which are 
the mobiles covered by sector antenna Sj but associated with 
a cell sector other than X; . 



Considering intracell interference, all of the mobiles within 
the sector transmit with power given by (TTZt . Substituting (fT2] i 
and (O into (J6j), we obtain the normalized received power of 
the intracell interferers 
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Although the number of mobiles in the reference cell sector 
must be known to compute the outage probability, the locations 
of these mobiles in the cell are irrelevant to the computation 
of the fli j of the intracell interferers. 

Considering intercell interference, the set Aj\Xj can be 
further partitioned into sets Aj D Xk, k ^ j, containing the 
mobiles covered by sector antenna Sj but associated with some 
other sector antenna Sk- For those mobiles in Aj DXk, power 
control implies that 

PdO^ /1Q f (\\S k - Xi\\) = P , X.eXkHA,, k^j. 



Substituting (O, (JT2J), and © into © yields 
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for Aj\Xj, which gives the mean intercell interference power 
at the reference sector antenna due to interference from mobile 
i of sector k — g(i). Assuming that imperfect power control 
is caused by numerous independent multiplicative effects, it 
may be approximately modeled by increasing the standard 
deviations of £jj, £ rj , and ^fc. 

B. Rate Control 

In addition to controlling the transmitted power, the rate 
Ri of each uplink needs to be selected. Due to the irregular 
network geometry, which results in cell sectors of variable 
areas and numbers of intracell mobiles, the amount of inter- 
ference received by a sector antenna can vary dramatically 
from one sector to another. With a fixed rate for each sector, 
the result is a highly variable outage probability. An alternative 
to using a fixed rate for the entire network is to adapt the rate 
of each uplink to satisfy an outage constraint or maximize the 
throughput of each uplink. Many networks adapt the modula- 
tion and coding to increase the throughput in accordance with 
the channel conditions. For example, in an 802. llg network, 
more than ten combinations of modulation and coding are 
available, including QPSK, 16-QAM, and 64-QAM. 

To illustrate the influence of rate on performance, consider 
the following example. The network has C = 50 base stations 
and M = 400 mobile stations placed in a circular network of 
radius r net = 2. The base-station exclusion zones have radius 
rbs = 0.25, while the mobile exclusion zones have radius 
r m = 0.01. The spreading factor is G = 16, and the chip factor 
is h = 2/3. Since M/C = G/2, the network is characterized 
as being half loaded. The SNR is V = 10 dB, the activity factor 
is pi = 1, the path-loss exponent is a = 3, and shadowing is 




Fig. 2. Outage probability of eight randomly selected uplinks (dashed lines) 
along with the average outage probability for the entire network (solid line). 
The results are for a half-loaded network (M/C = G/2), with distance- 
dependent fading and shadowing (cr s = 8 dB) and are shown as a function 
of the rate R. 

assumed with a s = 8 dB. A distance-dependent fading model 
is assumed, where m^j is set according to: 



if ||^-Xi||<rb,/2 

in,., <, 2 if r b5 /2 < \\Sj - Xi\ < r bs 
if \\Sj-Xi\ >r b5 



(16) 



The distance-dependent-fading model characterizes the situa- 
tion where a mobile close to the base station is in the line-of- 
sight, while mobiles farther away tend to be non-LOS. 

Fig. [2] shows the outage probability as a function of 
rate. Assuming the use of a capacity-approaching code, two- 
dimensional signaling over an AWGN channel, and Gaussian 
interference, the SINR threshold corresponding to rate R is 
j3 = 2^ — 1, The dashed lines in Fig. [2] were generated 
by selecting eight random uplinks and computing the outage 
probability for each using this threshold. Despite the use of 
power control, there is considerable variability in the outage 
probability. The outage probabilities {e;} were computed 
for all M uplinks in the system, and the average outage 
probability, 



1 M 



(17) 



is displayed as a solid line in the figure. 

Fig. [3] shows the throughput as a function of rate, where the 
throughput of the i th uplink is found as 



Ri(l 



(18) 



and represents the rate of successful transmissions. The param- 
eters are the same as those used to produce Fig. |2] and again 
the SINR threshold corresponding to rate R is /3 = 2 R — 1. The 
plot shows the throughput for the same eight uplinks whose 
outage was shown in Fig. as well as the average throughput 
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for the number of uplink users that are served. By taking into 
account the number of mobiles per unit area, the total data 
flow in a given area can be characterized by the transmission 
capacity, defined as lfl2l 



Fig. 3. Throughput of eight randomly selected uplinks (dashed lines) 
along with the average throughput for the entire network (solid line). System 
parameters are the same used to generate Fig. [2] 



1) Fixed-Rate Policies: Consider the case that all uplinks in 
the system must use the same rate; i.e., Ri = R, for all uplinks. 
On the one hand, the rate could be selected to maximize the 
average throughput. With respect to the example shown in Fig. 
12 this corresponds to selecting the R that maximizes the solid 
curve, which occurs at R — 1.81, However, at the rate that 
maximizes throughput, the corresponding outage probability 
could be unacceptably high. When R = 1.81 in the example, 
the corresponding average outage probability is E[e] = 0.37, 
which is too high for many applications. As an alternative to 
maximizing throughput, the rate R could be selected to satisfy 
an outage constraint £ so that E[e] < £. For instance, setting 
R = 0.84 in the example satisfies an average outage constraint 
£ = 0.1 with equality. To distinguish between the two fixed- 
rate policies, we call the first policy maximal-throughput fixed 
rate (MTFR) and the second policy outage-constrained fixed 
rate (OCFR). 

2) Variable-Rate Policies: If R is selected to satisfy an 
average outage constraint, the outage probability of the in- 
dividual uplinks will vary around this average. Furthermore, 
selecting R to maximize the average throughput does not 
generally maximize the throughput of the individual uplinks. 
These issues can be alleviated by selecting the rates Ri 
independently for the different uplinks. The rates could be 
selected to maximize the throughput of each uplink; i.e., Ri = 
arg max T.- L for each uplink, where the maximization is over all 
possible rates. We call this the maximal-throughput variable- 
rate (MTVR) policy. Alternatively, the selection could be made 
to require all uplinks to satisfy the outage constraint £; i.e., 
ti < C f° r a ll i< We call this the outage-constrained variable- 
rate (OCVR) policy. 

V. Performance Analysis 

A. Performance Metrics 

While the outage probability, throughput, and rate charac- 
terize the performance of a single uplink, they do not quantify 
the total data flow in the network because they do not account 
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(20) 



where A = M/irr^ is the density of transmissions in the 
network, and the units are bits per channel use per unit 
area. Transmission capacity can be interpreted as the spatial 
intensity of transmissions; i.e., the rate of successful data 
transmission per unit area. 

Performance metrics can be obtained using a Monte Carlo 
approach as follows. Draw a realization of the network by 
placing C base stations and M mobiles within the disk of 
radius r ne t according to the uniform clustering model with 
minimum base-station separation rb s and minimum mobile 
separation r m . Compute the path loss from each base station to 
each mobile, applying randomly generated shadowing factors. 
Determine the set of mobiles associated with each cell sector. 
Assuming that the number of mobiles served in a cell sector 
cannot exceed G, which is the number of orthogonal sequences 
available for the downlink, then set the rate of the last Mj — G 
mobiles in the cell sector (if there are that many) to zero (since 
they will be denied access). At each sector antenna, apply 
the power-control policy to determine the power it receives 
from each mobile that it serves. For each cell sector, compute 
the rates, outage probabilities, and throughputs according to 
the MTFR, OCFR, MTVR, or OCVR network policy. Then 
compute the transmission capacity. 

B. Simulation Parameters 

In the following subsections, the Monte Carlo method 
described in the previous subsection is used to characterize 
the uplink performance. In all cases considered, the network 
has C = 50 base stations placed in a circular network of 
radius r net = 2. Except for Subsection IV-EI which studies 
the influence of r^ s , the base-station exclusion zones are 
set to have radius rb s = 0.25. A variable number M of 
mobiles are placed within the network using exclusion zones 
of radius r m = 0.01. The SNR is T = 10 dB, and the 
activity factor is p. L = 1. Unless otherwise stated, the path-loss 
exponent is a = 3. The fading is the distance-dependent fading 
specified by ( TTol l. Both unshadowed and shadowed (<j s = 8 
dB) environments are considered. The chip factor is h = 2/3, 
and except for Subsection I V-DI which studies the influence of 
G, the spreading factor is G = 16. 

C. Policy Comparison 

Fig. H shows the average transmission capacity of the four 
network policies in distance-dependent fading, both with and 
without shadowing, as a function of the load M/C. For the 
OCFR policy, the average was obtained by optimizing over 
a rate that is common to all sectors. For the OCVR policy, 
the uplink rate Ri of each uplink was maximized subject to 
an outage constraint. While the transmission capacities of the 
MTFR and MTVR policies are potentially superior to those 
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Transmission capacity for the four network policies as function of 
M/C for distance-dependent fading and both shadowed (<r s = 8 
unshadowed cases, assuming that the rates are optimized. 



3.5- 

3 - 
2.5 - 



-▼— MTVR: M/C=G 

MTVR: M/C=G/2 
-•— MTFR: M/C=G 

■ • - MTFR: M/C=G/2 
-e— OCVR: M/C=G 

■ e - OCVR: M/C=G/2 
OCFR: M/C=G 

■ OCFR: M/C=G/2 




Fig. 5. Transmission capacity as function of spreading factor G for two 
values of system load, distance-dependent fading, and shadowing with cr s = 8 
dB. 

of the OCFR and OCVR policies, this advantage comes at 
the cost of a variable and high value of e, which is generally 
too large for most applications. The bottom pair of curves in 
Fig. [4] indicate that the OCVR policy has a higher average 
transmission capacity than the OCFR policy. 

D. Spreading Factor 

Fig. |5] shows the transmission capacity as a function of the 
spreading factor G (with h = 2/3) for the shadowed distance- 
dependent-fading channel. Two loads are shown for each of 
the four policies. An increase in G is beneficial for all policies, 
but the MTVR policy benefits the most. 

E. Base-Station Exclusion Zone 

Fig. [6] shows the transmission capacity for each of the four 
policies as a function of the base-station exclusion zone r^s 
for M/C — G/2 and two values of path-loss exponent a. 
The distance-dependent fading model is used, and shadowing 
is applied with cr s = 8 dB. The two policies that constrain the 
outage probability are more sensitive to the value of r^ s than 
the two policies that maximize throughput. 
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Fig. 6. Transmission capacity as a function of the base-station exclusion 
zone rt, s for four policies and two values of path-loss exponent a. 

VI. Conclusion 

A new analysis of DS-CDMA uplinks has been presented. 
This analysis is much more detailed and accurate than existing 
ones and facilitates the resolution of network design issues. 
In particular, it has been shown that once power control is 
established, the rate can be allocated according to a fixed-rate 
or variable-rate policy with the objective of either maximizing 
throughput or meeting an outage constraint. An advantage 
of variable-rate power control is that it allows an outage 
constraint to be enforced on every uplink, which is impossible 
when a fixed rate is used throughout the network. Another 
advantage is an increased transmission capacity. 
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